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Abstract 

Genetically controlled self-incompatibility systems represent links between 
genetic diversity and plant demography with the potential to directly impact on 
population dynamics. We use an individual-based spatial simulation to investi- 
gate the demographic and genetic consequences of different self-incompatibility 
systems for plants that vary in reproductive capacity and lifespan. The results 
support the idea that, in the absence of inbreeding effects, populations of self- 
incompatible species will often be smaller and less viable than self-compatible 
species, particularly for shorter-lived organisms or where potential fecundity is 
low. At high ovule production and low mortality, self-incompatible and self- 
compatible species are demographically similar, thus self-incompatibility does 
not automatically lead to reduced mate availability or population viability. 
Overall, sporophytic codominant self-incompatibility was more limiting than 
gametophytic or sporophytic dominant systems, which generally behaved in a 
similar fashion. Under a narrow range of conditions, the sporophytic dominant 
system maintained marginally greater mate availability owing to the production 
of S locus homozygotes. While self-incompatibility reduces population size and 
persistence for a broad range of conditions, the actual number of S alleles, 
beyond that required for reproduction, is important for only a subset of life 
histories. For these situations, results suggest that addition of new S alleles may 
result in significant demographic rescue. 



Introduction 

A central theme of evolutionary ecology is to understand 
how the interplay of genetic and demographic factors 
influences population dynamics, species persistence, and 
evolutionary processes. Quantifying the relative impor- 
tance of such factors also has direct application in many 
areas. For example, management of small plant popula- 
tions is a major concern in conservation biology with 
recent research efforts focusing on quantitative analysis of 
population viability. Major constraints considered to limit 
population performance in these situations fall into two 
broad classes. Demographic factors include asynchrony of 
flowering, pollinator limitation (Jennersten 1988; Ooster- 
meijer et al. 2003), demographic and environmental sto- 
chasticity (Holsinger 2000), and Allee effects (Stephens 
and Sutherland 1999). Genetic factors that have been 
implicated include mutational meltdown (Lynch and 
Gabriel 1990; Lynch et al. 1995), inbreeding depression 



(Oostermeijer 2000; Frankham 2005), and hybridization 
(Rhymer and Simberloff 1996; Field et al. 2009). 

It has been suggested that for those plant species with 
genetically controlled self-incompatibility (SI) systems 
(Richards 1997), reduced mate availability owing to loss of 
allelic richness (S alleles) at the incompatibility locus may 
also impose significant demographic constraints on small 
populations (Barrett and Kohn 1991; Byers and Meagher 
1992; Frankham et al. 2002; Castric and Vekemans 2004), 
especially when acting in concert with ecological factors 
such as pollinator limitation (Busch and Schoen 2008; 
Young et al. 2012). Such direct genetically mediated 
demographic constraints may represent one of the few 
clear links between genetic diversity and population viabil- 
ity of immediate importance for plant population dynam- 
ics and long-term persistence. As such, self-incompatibility 
represents an example of a genetic trait whose phenotypic 
expression directly influences an ecological process 
(Hughes et al. 2008). 
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Genetically controlled, homomorphic, self-incompati- 
bility occurs in approximately 50% of angiosperm families 
(Richards 1997; Igic et al. 2008). Allelic genealogy studies 
suggest that the two broad types of self-incompatibility 
systems, those exhibiting gametophytic control of mating 
phenotype, and those in which mating reactions are gov- 
erned by the sporophytic generation, have evolved inde- 
pendently and possibly several times (Matton et al. 1994). 
Presumably, this convergent evolution reflects the advan- 
tages that such mechanisms provide in preventing 
inbreeding and its associated negative fitness effects. 

Self-incompatibility systems function effectively in large 
populations where negative frequency-dependent selection 
acts to maintain large numbers of S alleles resulting in 
high mate availability. Emerson (1939) was the first to 
report the numbers of S alleles in a population, estimat- 
ing that there were 45 different alleles in a population of 
Oenothera organensis (which exhibits gametophytic SI) 
consisting of 500 individuals. Since then, self-incompati- 
bility loci have generally been shown to exhibit high 
within-population allelic richness (Richman and Kohn 
1996; Lawrence 2000; Edh et al. 2009). In small popula- 
tions, however, where founder events and genetic drift 
can erode allelic richness despite selection, self-incompati- 
bility systems may limit mate availability and, through 
this, fertilization success and fecundity (Busch and Schoen 
2008). This is more likely for sporophytic systems in 
which dominance relationships among alleles further 
reduce the effectiveness of frequency-dependent selection. 

Several empirical studies have invoked low S allele rich- 
ness as a possible cause of reduced seed set in small pop- 
ulations of self-incompatible plants: Hymenoxys acaulis 
var. glabra (DeMauro 1993), Helenium virginicum (Mess- 
more and Knox 1997), Cochlearia bavarica (Fischer et al. 
2003), Ranunculus reptans (Willi et al. 2005), and Linnaea 
borealis (Scobie and Wilcock 2009). Recently, Young and 
Pickup (2010) found positive relationships between 
population size, S allele richness, and mate availability in 
populations of the endangered daisy Rutidosis leptorrhyn- 
choides that has a sporophytic self-incompatibility system. 
In this case, populations with few S alleles also exhibited 
elevated interplant variance in male and female fitness 
components, as expected under mate limitation (Byers 
and Meagher 1992), and reduced long-term viability 
as assessed by matrix projection models (Young et al. 
2000). 

Both analytical and computer simulation models have 
been used to examine how population size and gene flow 
interact to determine S allele diversity (Imrie et al. 1972; 
Schierup 1998; Castric and Vekemans 2004; Kirchner 
et al. 2006; Cartwright 2009; Levin et al. 2009) and how 
this affects mate availability and reproductive success in 
subdivided populations (Kirchner et al. 2006; Wagenius 



et al. 2007; Levin et al. 2009). Vekemans et al. (1998) 
used both deterministic and stochastic simulation 
approaches to show that gametophytic systems are less 
mate limiting than sporophytic systems and that domi- 
nance relationships among S alleles in sporophytic sys- 
tems, and multiple pollination events strongly affect mate 
availability. Metapopulation models also suggest that 
effective among-population migration rates are likely to 
be higher for dominant than recessive alleles at SI loci 
(Schierup et al. 2000). 

The combined results of these empirical and simulation 
investigations suggest that S allele number, frequency, and 
dominance relationships may have significant effects on 
the demography and genetic structure of plant popula- 
tions through their influence on mating opportunities, 
and should be considered as important determinants of 
population viability for self-incompatible species. Some 
theoretical studies have begun to illustrate how demo- 
graphic and genetic processes can interact to affect popu- 
lation viability (Kirchner et al. 2006; Wagenius et al. 
2007; Levin et al. 2009). Here, we extend these analyses 
to explicitly investigate the influence of life-history and 
reproductive traits on individual fitness with the aim of 
elucidating broad predictive relationships that may drive 
population performance. For example, the relative impact 
of different mating systems is likely to vary with plant 
life-history features (e.g., birth and death rates) and this 
has so far not been considered. 

Here, we develop a spatially explicit individual-based 
simulation approach which extends previous empirical 
and modeling studies to simultaneously examine genetic 
and demographic effects of both gametophytic and sporo- 
phytic SI on population viability. The important point of 
difference between these two SI systems is that in the 
gametophytic system, the pollen mating phenotype is 
controlled by the haploid pollen genotype, while in the 
sporophytic system, the pollen mating phenotype is deter- 
mined by the diploid male genotype. Sporophytic SI 
therefore generally results in more limited mate availabil- 
ity although this may be alleviated by the effects of domi- 
nance relationships among S alleles. In addition, we assess 
how the relative impact of these systems is influenced by 
variation in major life-history parameters and the diver- 
sity of alleles at the SI locus. This integrated approach has 
several important advantages. First, we can identify sets of 
conditions where demographic or genetic factors are more 
limiting, which has practical implications for conservation 
management. Second, model predictions can be explicitly 
linked to empirically measureable population parameters 
(e.g., allelic richness, levels of correlated paternity, vari- 
ance in seed set). Finally, it may be possible to identify 
suites of life-history traits that make species more demo- 
graphically sensitive to changes in S allele richness. 
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Computer Simulations 
Biological assumptions 

To investigate population viability in relation to self- 
incompatibility and life history, the model allows plants to 
vary in reproductive capacity and lifespan (i.e., annual to 
long-lived perennials). Seed formation following a mating 
event is determined by a single SI locus with variable num- 
bers of alleles. In the results presented here, this was varied 
from three (the minimum needed for there to be at least 
the possibility of successful mating events to occur) to 50, 
which approaches the maximum level of diversity seen at 
the population level in many species (Lawrence 2000; 
Castric and Vekemans 2004). The model also assumes five 
unlinked variable neutral loci, each with five alleles. In 
addition, the model assumes diploid genetics and that indi- 
viduals can act both as pollen donors and as females. 

The primary aim of this study was to gain understanding 
of how these major life-history features jointly influence 
variation in ecological factors associated with individual fit- 
ness (e.g., mean and variance in numbers of compatible 
mates and pollen donors, per capita seed production), as 
well as to examine how shifts in the mating system might 
interact with potential fecundity and death rate to impact 
on frequently measured genetic parameters (e.g., allelic 
richness at the SI locus and neutral loci, observed [H 0 ] and 
expected [H E ] heterozygosity, and the inbreeding coeffi- 
cient [-FisD- We therefore use an individual-based approach 
to characterize variables that are frequently measured in 
field studies and investigate how these relate to population 
viability and persistence. The model assumes no fitness 
consequences of selfing (i.e., when using a self-compatible 
mating model as a baseline comparison), as the goal is not 
to examine the conditions under which self-incompatibility 
evolves, but rather to focus on the implications of SI for 
plant population viability. 

General model structure and initial 
conditions 

The simulations (written in C++ using the GNU Scientific 
Library) consist of a two-dimensional (100 x 100) grid, 
with absorbing boundaries, as this scenario best mimics 
natural populations of finite size, beyond which propa- 
gules are essentially lost from the system. The fraction of 
sites suitable for occupancy can be varied - these are 
assigned randomly at the beginning of each run (this 
value was fixed at 0.25 for the results reported here). At 
the start of each run, genotypes were randomly assigned 
to each of 200 individuals (with the constraint that indi- 
viduals must be heterozygous at the SI locus). These were 
randomly placed in a subset of the available spatial 



locations. In addition to spatial location and genotype, 
the simulation also tracks individual ages. 

In each time interval, for each reproductive adult (the 
minimum age of reproduction was fixed at 1, thus all indi- 
viduals were considered to be adults), the total number of 
possible pollen donors (including both compatible and 
incompatible) within pollination distance is determined. 
The spatial scales of both pollen and seed movement can 
be varied in the simulation, however, for the purposes of 
this study, these were fixed at 10 and 5, respectively (e.g., 
for each adult, potential pollen donors were chosen from 
the pool of individuals within 10 units or less). 

To calculate per capita ovule production, we assumed a 
linear relationship with individual age (A), such that 
M = the age of first reproduction, L = the expected life- 
span {lid, where d = death rate), and b 0 = maximum per 
capita ovule production. The number of ovules produced 
by an individual (Ov) is thus {b 0 d} A when A is less than 
or equal to I, otherwise Ov = b 0 (this is reasonable for 
most plants, particularly long-lived perennials, where 
reproductive output does not increase indefinitely with 
size). To ensure that the number of ovules produced is at 
least 1 at the minimum reproductive age, b 0 was con- 
strained to be greater than 1/dM. For each ovule pro- 
duced by an individual, a random pollen donor was 
chosen from the pool of possible mates (haploid ovule 
and pollen genotypes were randomly constructed from 
the diploid parents), and then, compatibility was deter- 
mined based on whether the mating system was self- 
compatible (i.e., full compatibility across all individuals), 
gametophytic SI, sporophytic codominant, or sporophytic 
dominant (dominance was assumed to be paternally con- 
trolled as this is the most commonly seen form; Richards 
1997). In the case that the ovule and the chosen mate are 
incompatible, then the ovule is aborted and another ovule 
is chosen for mating. For each plant, seeds generated 
from successful fertilizations were then dispersed ran- 
domly within the seed dispersal range. This approach 
simulates pollination by multiple pollen donors at the 
maternal level, which has been shown to be important in 
generating biologically realistic fertilization probabilities 
(Vekemans et al. 1998). Seeds dispersing outside of the 
grid or to cells already occupied by adults were assumed 
to have been lost. Once all individuals had mated and 
reproduced, and seeds dispersed to empty cells, adult 
death occurred with a fixed probability. Following this, in 
each cell where one or more seeds had landed, a single 
individual was randomly chosen to germinate. 

Simulation experiments 

To explore the demographic and genetic consequences of 
variation in self-incompatibility system, S allele richness 
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and life history, the following parameters were varied sys- 
tematically in the simulation model: mating system (as 
noted above), the number of S alleles (3, 5, 10, 25, 50), 
death rate (0.05, 0.1, 0.2, 0.3, 0.5, 1.0), and per capita 
ovule production rate (10, 25, 50, 75, 100). Note that 
even very high ovule production values translated into 
realistic numbers for per capita production of surviving 
offspring. The situations we investigated simulated life 
histories ranging from biennials to long-lived perennials 
and covered a wide range of reproductive values. 

Each simulation ran for a maximum of 500 generations - 
for runs in which the population persisted for the full 
duration, both demographic and genetic patterns gener- 
ally stabilized within the first 100-200 generations. For 
each set of the parameter combinations outlined above, 
100 random runs were conducted, and results averaged 
across runs for the following demographic variables at 
each generation: total number of reproductive plants, 
proportion of reproductive individuals within pollen-flow 
range that were of compatible genotypes at the SI locus 
(mate availability), number of seeds/plant and its vari- 
ance. Mate availability was calculated as the average num- 
ber of compatible individuals. Similarly, information was 
gathered at each generation for several genetic variables: 
number of alleles per neutral locus (across loci and indi- 
viduals), number of alleles at the SI locus, observed and 
expected heterozygosity, and the fixation index (Fi S ). In 
addition, for each set of parameter combinations, the per- 
centage of runs that persisted for the full 500 generations 
was recorded, as well as average population persistence 
time (number of generations). 

Results 

Population viability 

Negative effects on population persistence were generally 
more dramatic with increasing death rates and lower 
birth rates. Regardless of which self-incompatibility sys- 
tem was assumed, or the level of allelic richness at the SI 
locus, when death rates were 0.3 or greater, populations 
uniformly went extinct within a few generations (<50). 
For lower death rates (d < 0.2), there was significant vari- 
ation among SI systems with respect to persistence and 
population size, and clear effects of genetic diversity at 
the SI locus (Figs. 1, 2). Regardless of death rate, ovule 
production or S allele richness, the sporophytic codomi- 
nant system always resulted in the lowest population per- 
sistence and size. The gametophytic and sporophytic 
dominant systems generally behaved in a similar fashion, 
with persistence times and population sizes intermediate 
to the sporophytic codominant and self-compatible 
situations. 



P. H. Thrall ef a/. 

Time to extinction 

For very low death rates (d = 0.05), population persis- 
tence was only significantly reduced relative to the self- 
compatible case when ovule production was also very low 
(b Q < 25). S allele richness had no effect on population 
persistence for either the gametophytic or sporophytic 
dominant systems. However, for the sporophytic codomi- 
nant system, populations initiated with only five S alleles 
had an increased probability of going extinct for substan- 
tially higher rates of ovule production (b Q < 50; Fig. 1J). 

For somewhat higher death rates (d = 0.1-0.2), limita- 
tions on population persistence imposed by self-incom- 
patibility were manifested over a broader range of ovule 
production rates and levels of S allele richness. At the 
lowest level of ovule production (b a = 10), these death 
rates also led to reduced population persistence time for 
the self-compatible case. Interactions between S richness, 
ovule production, and mating system also became appar- 
ent. These were most obvious for the sporophytic codom- 
inant system where the effect of reduced S allele richness 
was strongly dependent on ovule production. For exam- 
ple, when d = 0.1 and the ovule production rate was 50, 
reducing S allele richness from 25 to 5 decreased persis- 
tence time by 80% (i.e., from 500 to 100 generations; 
compare Fig. IE and K). For the same death rate, but 
with 75 ovules produced per plant, then the effect of 
reducing S allele richness was only a 50% drop in persis- 
tence, and no reduction was observed when the ovule 
production rate was 100. This interaction between ovule 
production and S allele richness was evident to a lesser 
extent for the gametophytic and sporophytic dominant 
systems. 

Further increases in death rate (d = 0.2) resulted in a 
general and dramatic reduction in population persistence 
for all self-incompatibility systems, with long-term popu- 
lation viability only being achieved for the highest rates of 
ovule production and then only for the gametophytic and 
sporophytic dominant systems. These two SI systems were 
also strongly influenced by S allele richness but, in con- 
trast to the case when d = 0.1, effects were most evident 
at higher ovule production rates (b Q = 50-75; Fig. 1). 
This was not the case for the sporophytic codominant 
system, where persistence was uniformly low, and extinc- 
tion rates were little affected by either S allele richness or 
ovule production. 

Population size 

Although, in the simulation results presented here, a fixed 
2500 of the 10,000 sites represented suitable habitat, in 
no case did the percentage of sites occupied exceed 50% 
of the available habitat. This indicates that constraints on 
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Death = 0.05 (B) Death = 0.10 (C) Death = 0.20 



Figure 1. Population persistence under 
different breeding systems as a function of 
death rate, ovule production rate, and number 
of alleles at the self-incompatibility (SI) locus. 
Moving from left to right for the three 
columns of graphs, death rate was fixed at 
d= 0.05, 0.10, and 0.20, respectively. Moving 
vertically down each column, the number of 5 
alleles was fixed at 50, 25, 10, and 5, 
respectively. For each graph, ovule production 
rate was varied (b D = 10, 25, 50, 70, 100). 
Each graph shows persistence for the 
gametophytic (solid orange lines), sporophytic 
dominant (solid blue lines), and sporophytic 
codominant (solid purple lines) SI systems. 
Note that gametophytic and sporophytic 
dominant SI systems show very similar 
persistence for d = 0.05 and 0.1, and in these 
cases, the lines overlap. Results for the self- 
compatible case are plotted on all graphs for 
comparison (dashed lines). 




50 100 



Ovules plant 



population growth were demographic and genetic, rather 
than being due to the simulation structure itself (i.e., 
there were always available sites for plants to colonize). 
Trends in population size after 500 generations paralleled 
the patterns seen with respect to persistence, although the 
effects of decreases in ovule production were more 
pronounced, even for the self-compatible case (Fig. 2). 

Irrespective of death rate, for the self-compatible case, 
there was a clear asymptotic relationship between ovule 
production and population size, such that beyond a cer- 
tain level, further increases in ovule production did not 
result in a higher percentage of the available sites being 
occupied by adults. Moreover, this asymptote itself 
decreased with increasing death rates. For low death rates 
(d = 0.05, 0.1), this pattern was also evident for the three 
SI mating systems, but significant reductions in popula- 
tion size relative to the self-compatible case were observed 
when S allele richness was 10 or less and ovule production 



lower than 75, particularly for the codominant sporo- 
phytic case. As with persistence time, the interacting 
effects of S allele richness and ovule production were most 
evident when d = 0.1. When d = 0.2, the magnitude of 
reductions in population size relative to the self-compati- 
ble case was great enough that the asymptote would only 
be reached for ovule production rates greater than 100, 
regardless of which SI mating system was assumed. 

Taken together, these results indicate that at very low 
death rates (d = 0.05), population sizes are generally large 
enough that neither demographic nor genetic factors 
impact significantly on persistence, except when both 
ovule production rates and S allele richness are very low. 
In contrast, at death rates of 0.2 or greater, demographic 
limitations reduce population size to the extent that even 
with high ovule production rates, large numbers of S 
alleles are necessary for long-term persistence. However, 
even in this case, population sizes remain small. It is only 
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(A) Death = 0.05 (B) Death = 0.10 (C) Death = 0.20 




Figure 2. Population size after 500 
generations under different breeding systems 
as a function of death rate, ovule production 
rate, and number of alleles at the self- 
incompatibility (SI) locus. Moving from left to 
right for the three columns of graphs, death 
rate was fixed at d = 0.05, 0.10, and 0.20, 
respectively. Moving vertically down each 
column, the number of S alleles was fixed at 
50, 25, 10, and 5, respectively. For each 
graph, ovule production rate was varied 
(b 0 = 10, 25, 50, 70, 100). Each graph shows 
population size for the gametophytic (solid 
orange lines), sporophytic dominant (solid blue 
lines), and sporophytic codominant (solid 
purple lines) SI systems. Results for the self- 
compatible case are plotted on all graphs for 
comparison (dashed lines). Note that 
gametophytic and sporophytic dominant SI 
systems show very similar population sizes for 
d = 0.05 and 0.1, and in these cases, the lines 
overlap. Population size is expressed as the 



Ovules plant 1 proportion of suitable sites occupied. 



at intermediate death rates (d = 0.1) where population 
sizes and persistence can be strongly influenced by both 
demographic and genetic factors. 

Demographic and genetic dynamics 

Changes in persistence and equilibrium population size 
are the most visible consequences of the impacts of shifts 
in major life-history parameters (e.g., birth and death 
rates) on a range of underlying demographic (mate avail- 
ability, seed set, male/female fitness) and genetic (loss of 
S alleles, H 0 , F IS ) factors. In self-incompatible mating sys- 
tems, these variables can interact to severely reduce popu- 
lation sizes, thereby resulting in decreased long-term 
persistence through stochastic extinction events. The 
results presented above indicate that demographic and 
genetic effects are both important for intermediate death 
rates and low to intermediate levels of fecundity and S 



allele richness. To examine in detail how these interact to 
impact on viability, in this section, we examine popula- 
tion dynamics for situations where d = 0.1, b 0 < 75, and 
the number of S alleles <50. We focus on comparisons 
between the gametophytic and sporophytic codominant 
situations as these bracket the extremes in terms of how 
restrictive SI systems can be, because, with few exceptions 
(see Gametophytic vs. Sporophytic Dominant Self-Incom- 
patibility below), the sporophytic dominant and gameto- 
phytic cases were generally similar. 

Mate availability and F, s 

When ovule production was very low (b Q = 10), S alleles 
were rapidly lost from populations regardless of breeding 
system or initial S allele richness. In these situations, this 
reduction resulted in a severe restriction in local mate avail- 
ability relative to the self-compatible case, very negative 



678 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



P. H. Thrall ef al. 



Effects of Self-incompatibility on Viability 



inbreeding coefficients (F IS ) indicating strong disassortative 
mating (Fig. 3A and B), and extinction always occurred 
within 150 generations (Fig. 1). 

Relatively small increases in ovule production (b Q = 25) 
were sufficient to ensure population growth and long-term 
viability of self-compatible populations but not for either 
the gametophytic and sporophytic codominant SI systems, 
for which population sizes still declined substantially from 
initial values (Fig. 3C and D). However, the influence of S 
allele richness on mate availability and F ls was now depen- 
dent on breeding system. For the gametophytic system, 
increasing S allele diversity led to reduced disassortative 
mating and greater retention of mate availability (Fig. 3C). 
In contrast, the sporophytic codominant system continued 
to be severely mate-limited, regardless of the number of S 
alleles with which the simulation was initiated. Further 
increases in ovule production (b Q = 50) resulted in positive 
population growth for the gametophytic system regardless 
of S allele richness - mate availability never declined and 
F IS was always either zero or slightly positive. However, S 
allele richness had dramatic effects on these parameters for 
outcomes in the sporophytic codominant case (Fig. 3F), 
ranging from extreme mate limitation and extinction 
(S = 10) to random mating and demographic increase 
similar to the self-compatible situation (S = 25). 

Seed set and variance in fitness 

Changes in mate availability associated with population 
size jointly depend on breeding system and S allele diver- 
sity. In turn, the former exert their influence on popula- 
tion growth and persistence primarily through their 
effects on individual reproductive performance. This is 
most easily observed by examining the mean and variance 
of seed set. 

In situations where low numbers of S alleles or low 
ovule production rates led to severe mate limitation (e.g., 
Fig. 4A, B, and D), then decline in population size was 
correlated with an initial increase in both the mean and 
variance of seed set. This phase of the dynamics was fol- 
lowed by a rapid decline in both variables as disassortative 
mating reduced mate availability to very low levels. This 
partly reflects the early development of an age-structured 
population (ovule production rates were assumed to be a 
positive function of plant age/size) and subsequently the 
buildup of fine-scale spatial genetic structure. Such struc- 
ture initially provides increased opportunities for mating 
as local plant density increases, but as S alleles are lost and 
relatedness increases, the local pool of genetically compati- 
ble pollen donors ultimately becomes restricted. 

The point at which the positive effect of increasing S 
alleles on mate availability is sufficient to maintain seed 
set and prevent extinction (Fig. 4C, E, and F) depends on 



breeding system and the maximum rate of ovule produc- 
tion. Higher rates of ovule production and gametophytic 
systems require fewer S alleles for population persistence 
than lower ovule production rates and sporophytic 
codominant systems. Just above these thresholds (Fig. 4C), 
the maintenance of mate availability in small populations 
is not reflected in a qualitative difference in the general 
nature of the dynamical relationship between population 
size and seed set. Instead, persistence of small populations 
when S = 25 is due to a continued low level of seed set. 
As ovule production increases further (and demographic 
constraints become less limiting), both qualitative and 
quantitative effects of S allele diversity on the relationship 
between population size and mean and variance in fecun- 
dity become apparent (Fig. 4F). In this case, S = 25 popu- 
lations grow, while S = 10 populations persist at reduced 
size and S = 5 populations consistently go extinct. 

Gametophytic versus sporophytic dominant 
self-incompatibility 

Previous theoretical studies (Vekemans et al. 1998) have 
suggested that sporophytic self-incompatibility systems 
should generally be more limiting in terms of mate avail- 
ability than gametophytic systems due to the increased 
likelihood of matching S alleles between mates in sporo- 
phytic systems where the pollen mating phenotype is 
determined by the diploid male SI genotype. Our results 
indicate that under most conditions, these systems per- 
form in a very similar fashion (although if there were dif- 
ferences, it was usually the sporophytic dominant system 
that was more limiting). However, at very low levels of S 
allele diversity (where population sizes are also small), our 
results demonstrated that mate availability was actually 
higher for the sporophytic dominant system than for the 
gametophytic system (Fig. 5). The most likely explanation 
for this is that dominance relationships among alleles in 
the sporophytic system permit the buildup of SI locus 
homozygotes. This is not possible with either gameto- 
phytic or sporophytic codominant self-incompatibility 
systems. 

Discussion 

Effects of self-incompatibility on population 
viability 

Theoretical and simulation studies have suggested that 
self-incompatible breeding systems have significant poten- 
tial to negatively impact long-term population fecundity 
and viability, particularly during population crashes 
and subsequent bottlenecks, or in the colonization phase 
of population growth (e.g., Byers and Meagher 1992; 
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Figure 3. Dynamical relationship between population size, mate availability, and the inbreeding coefficient (F| S ) for different levels of S allele 
richness. The left-hand column of graphs shows results for the gametophytic self-incompatibility system, and the right-hand column of graphs 
shows results for the sporophytic codominant system. Ovule production rates were fixed at 10, 25, and 50 for the top, middle, and bottom rows 
of graphs, respectively. Changes in population size equate directly with dynamics through time, with arrows indicating the temporal direction of 
changes from the starting point (to). Trajectories represented by filled green circles are for S = 50, while filled blue and red circles represent 
dynamics for S = 25 and S = 10, respectively. Open black circles show dynamics under the self-compatible scenario. For all results shown, the 
death rate (of) was fixed at 0.1 . 
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Figure 4. Dynamical relationship between population size, average seed set, and variance in seed set for females for different levels of S allele 
richness. The left-hand column of graphs shows results for the gametophytic self-incompatibility system, and the right-hand column of graphs 
shows results for the sporophytic codominant system. Ovule production rates were fixed at 10, 25, and 50 for the top, middle, and bottom rows 
of graphs, respectively. Changes in population size equate directly with dynamics through time, with arrows indicating the temporal direction of 
changes from the starting point (f 0 ). Trajectories represented by filled green circles are for S = 50, while filled blue and red circles represent 
dynamics for 5= 25 and 5= 10, respectively. Open black circles show dynamics under the self-compatible scenario. For all results shown, the 
death rate (d) was fixed at 0.1. 
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Figure 5. Comparison of relationship between 
population size and mate availability for 
populations with gametophytic (solid orange 
line) and sporophytic (solid blue line) dominant 
self-incompatibility under different levels of 5 
allele richness: 5 = 5 (open circles); S = 10 
(open squares); 5 = 25 (diamonds). Results 
show that when the number of 5 alleles is ten 
or greater, gametophytic systems have higher 
mate availability, but when there are fewer 
than five S alleles, sporophytic incompatibility 
is less limiting. This is probably owing to the 
generation of 5 locus homozygotes in the 
sporophytic system that cannot arise under 
gametophytic control. The self-compatible 
situation is shown with open black circles and 
dashed lines. 



Vekemans et al. 1998; Wagenius et al. 2007; Hoebee et al. 
2008). In these situations, impacts are likely to be primar- 
ily through severe reductions in mate availability and 
fertilization success owing to reduced numbers of S 
alleles, as has been demonstrated in a number of single 
species case studies (Fischer et al. 2003; Willi et al. 2005; 
Glemin et al. 2008; Scobie and Wilcock 2009; Leducq 
et al. 2010; Young and Pickup 2010). 

Results from deterministic computer simulations of 
genetic dynamics indicate that in small populations, 
decreases in S allele diversity lead to reductions in average 
mate availability but an increase in the variance of avail- 
able mates (Byers and Meagher 1992; Vekemans et al. 
1998). Because these studies were focused explicitly on 
genetic rather than demographic processes, it is difficult 
to know whether these effects are likely to be significant 
relative to other population-level processes in natural sys- 
tems. However, recent studies suggest that both demo- 
graphic and genetic constraints operate in subdivided 
populations (Kirchner et al. 2006; Wagenius et al. 2007; 
Levin et al. 2009). The model we present explicitly inte- 
grates numerical and genetic dynamics in a biologically 
realistic spatial context, allowing us to directly evaluate 
the potential for self-incompatibility and life-history fea- 
tures to jointly impact on population viability relative to 
other demographic processes. Our results generally sup- 
port the idea that, in the absence of inbreeding effects, 
populations of self-incompatible species will often be 
smaller and less viable than self-compatible species, par- 
ticularly for shorter-lived organisms or where birth rates 
are low, although it must be remembered that the lack of 
inbreeding depression in the self compatible (SC) case 
may inflate this difference. Thus, one major conclusion is 
that longer-lived organisms such as trees can withstand neg- 
ative genetic effects on reproduction better than short-lived 



species. However, for some ovule production and death 
rates, SI systems performed in a similar manner demo- 
graphically to the self-compatible system, and thus, self- 
incompatibility does not automatically lead to reductions 
in mate availability. 

It has been hypothesized that sporophytic self-incom- 
patible systems should be more mate-limited than game- 
tophytic systems (e.g., Vekemans et al. 1998; Castric and 
Vekemans 2004), and some simulation studies have 
shown this (Levin et al. 2009). As expected, our results 
showed that, under most demographically limiting condi- 
tions, the sporophytic codominant system performed 
worse than either the sporophytic dominant or gameto- 
phytic systems - in the former case, more severe restric- 
tions on mate availability resulted in smaller populations 
that persisted for shorter periods of time (and in fact the 
available data indicate that fully codominant systems are 
rare in nature; Richards 1997). However, there were rela- 
tively few differences between the sporophytic dominant 
and gametophytic systems, although surprisingly under a 
limited set of conditions (e.g., when S allele richness was 
very low, S < 5), the sporophytic dominant system main- 
tained marginally greater mate availability due to the pro- 
duction of SI locus homozygotes, leading to small 
increases in population size and viability over the gameto- 
phytic system. Again in this comparison, the lack of 
inbreeding depression in the model may inflate the 
observed difference as it means that there is no fitness 
cost associated with highly homozygous (inbred) individ- 
uals while they do gain a mating advantage. 

5 allele richness and population viability 

Given that under some conditions, self-incompatibility 
can have significant viability consequences for plant 
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populations, it is important to ask when S allele richness 
(and the erosion of genetic variation at the SI locus) is of 
concern. This issue has been raised previously by several 
authors (Les et al. 1991; Young et al. 2000; Frankham 
et al. 2002). While possessing some form of self-incom- 
patibility clearly reduces population size and persistence 
for a broad range of conditions, the actual number of S 
alleles is important for a more limited set of life histories. 
For example, when death rates are high (and populations 
are severely demographically limited), then adding new S 
alleles has little or no effect on viability. This is also the 
case when demographic constraints are absent (low death 
rates and intermediate to high birth rates - for these 
parameter ranges, no genetic limitations on population 
size were apparent). It is only for intermediate situations, 
where differences between SI systems also become more 
apparent, that S allele richness further impacts on viability 
(and even then, this is generally only seen when the num- 
ber of alleles is relatively low). 

The available empirical evidence indicates that most 
populations of SI species maintain quite high numbers of 
S alleles (Richards 1997; Lawrence 2000; Castric and 
Vekemans 2004). Nevertheless, our results also indicate 
that under at least some conditions, increasing population 
viability through addition of new S alleles may be the most 
effective approach to population management (Fig. 6) and 
there is some empirical support for this strategy. An 
extreme example is the threatened daisy, Hymenoxys acau- 
lis var. glabra, which was shown to lack cross-compatible 



mating types in remaining populations; the recovery pro- 
gram explicitly focused on establishing new populations 
with a diversity of S allele types (DeMauro 1993). The 
potential for genetic rescue to increase fecundity through 
the addition of new S alleles has also been demonstrated 
experimentally for several other short-lived perennial 
herbs (Ranunculus reptans, Willi and Fischer 2005; Rutido- 
sis leptorrhynchoides, Pickup and Young 2008; Brassica 
insularis, Glemin et al. 2008). 

The effectiveness of such genetic rescue strategies partly 
depends on what type of self-incompatibility is present in 
a particular species. Based on theoretical considerations, it 
has been suggested that the efficiency of frequency-depen- 
dent processes in maintaining rare S alleles is likely to 
vary for the different SI systems (Vekemans et al. 1998; 
Levin et al. 2009). Our results (data not presented) 
showed that alleles were lost consistently slower in the 
gametophytic system than in either of the sporophytic 
systems (although overall the differences between sporo- 
phytic dominant and gametophytic SI were relatively 
small). The rate of loss (and the endpoint) also was 
dependent on the interaction between population size, 
mate availability, and seed set. 

Interactions between self-incompatibility, 
demography, and genetics 

Overall, our simulation results showed that, dynamically, 
population declines (as a consequence of high mortality 




Population size 



Figure 6. The influence of increasing the number of S alleles on population size in a species with sporophytic codominant self-incompatibility. 
When ovule production is low, increasing the number of S alleles from 5 to 10 has a strong positive influence on population size, indicating that 
genetic rescue is an effective management strategy (trajectory a - orange dashed arrow). Beyond this, however, increasing numbers of 5 alleles 
has much less effect (trajectory b - red dashed arrow), suggesting that demographic management will be more useful to promote population 
viability (trajectory c - blue dashed arrow). In this case, an effective increase in ovule production to d 0 = 100, achieved in practice through the 
promotion of pollination or germination, is approximately 25% more effective in promoting population recovery. 
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and/or low ovule production rates) are accompanied by a 
cascade of demographic and genetic effects. Some vari- 
ables respond only very slowly to decreases in population 
size (e.g., F 1S ) with little apparent change over large size 
ranges. Generally, F IS values were highest for intermediate 
sized populations, decreasing slightly for larger popula- 
tions (but always positive in these cases). As population 
sizes became very small f IS values often oscillated dramat- 
ically (sometimes between negative and positive values) 
before population extinction occurred, however, the gen- 
eral trend was toward negative values reflecting strong 
mate limitation and disassortative mating. The generation 
of significant negative values at the neutral loci as a result 
of disassortative mating suggests that physical linkage is 
not necessary for the S locus to exert considerable influ- 
ence on the genotypic makeup of individuals. This 
extends previous results showing that loci linked to the S 
locus are not only less amenable to purging of deleterious 
alleles by selection (Glemin et al. 2001), but are also pro- 
tected from the effects of random drift (Schierup et al. 
2000; Ruggiero et al. 2008). Note that under gameto- 
phytic self-incompatibility, populations could become 
much smaller before going extinct than with sporophytic 
self-incompatibility, particularly with respect to the 
codominant case. In the latter, mate availability became 
limiting far sooner. 

For other parameters (e.g., reduction of mate availabil- 
ity, variance in seed set), responses to changes in popula- 
tion size were much more predictable (i.e., we did not 
observe oscillations as for F IS ). Interestingly, both average 
mate availability and variance in seed set were lowest for 
very small or very large populations. When populations 
are large, then mate availability is not limiting and all 
individuals can set seed, while the extreme limitation at 
small sizes results in generally low fecundity. Initial 
increases in seed set, even though population sizes are 
declining (Fig. 4), are primarily due to increases in local 
plant density through recruitment into empty sites and 
the development of age structure. However, as spatial 
genetic structure builds up further, local clusters of indi- 
viduals become increasingly related (data not shown), 
particularly as S alleles are lost, and mate availability 
decreases again. Thus, not only are life-history character- 
istics important, but pollination and dispersal parameters 
are critical because of their impact on the buildup of local 
spatial genetic structure. Wells and Young (2002) present 
supporting evidence for this effect in the self-incompatible 
daisy, Rutidosis leptorrhynchoides. 

Results from empirical studies 

Several empirical studies have examined aspects of the 
population biology of self-incompatible species, most 



commonly limits placed upon reproductive performance 
by low genetic mate availability. For example, Les et al. 
(1991) investigated the reproductive biology and genetic 
structure of a rare daisy, Aster furcatus, which has sporo- 
phytic self-incompatibility and primarily persists in small 
populations. Their results showed that the inbreeding 
coefficient (Fi S ) was generally negative in these situations, 
indicating an excess of heterozygotes. Our simulation 
results suggest that this observation likely reflects strong 
disassortative mating under severe mate limitation when 
S allele numbers are low. Based on this result, Les et al. 
(1991) suggested that severe bottlenecks (e.g., through 
habitat fragmentation) leading to the loss of incompatibil- 
ity types may be a major threat to long-term persistence 
of species such as A. furcatus. Similar negative F 1S values 
have also been reported for several highly outcrossed 
tree species (Hall et al. 1994; Dayanandan et al. 1999; 
Degen et al. 1999; Lexer et al. 2000), indicating that this 
may be a common outcome for small populations of SI 
plants. 

Detailed quantitative analyses of mating structure, demog- 
raphy, and genetic variation have been carried out in multi- 
ple populations of several rare plant species with different 
life histories: a woody shrub (Grevillea iaspicula (Hoebee and 
Young 2001)) exhibiting gametophytic self- incompatibility 
and two herbaceous perennials showing sporophytic domi- 
nant self-incompatibility (Rutidosis leptorrhynchoides (Young 
et al. 1999, 2000); Arnica montana (Luijten et al. 2000)). 
Results from these studies provide strong support for the 
dynamical patterns seen in our model. For example, in both 
Arnica and Rutidosis, negative _F IS values were consistently 
observed in populations of <50 plants (Young et al. 1999), 
which is within the size range at which our model predicts 
this should occur. Young et al. (1999) suggested that nega- 
tive inbreeding coefficients might be due to an absence of 
homozygotes for low-frequency alleles. However, our analy- 
ses indicate that these patterns are much more likely to be 
due to severe mate limitation caused by the loss of S alleles. 
Further support for this interpretation comes from the 
observed high levels of correlated paternity and the strong 
differentiation between pollen and adult allele pools in small 
populations of Rutidosis (Young and Brown 1999; Young 
and Pickup 2010). Interestingly, the data presented in Young 
et al. (1999) also suggest that decreases in population size 
may have very little impact on F IS until a threshold is 
reached, below which values decline very sharply. This 
pattern closely matches dynamical predictions for _F IS derived 
from our model (see Fig. 3). 

Grevillea iaspicula is one of the few other self-incom- 
patible species whose demography, mating system, and 
genetic structure have been intensively studied. This spe- 
cies persists in a few small remnant sites with population 
sizes generally less than 50 flowering plants. Simulation 
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studies using the model described in this study again 
indicate close agreement with empirical observations of 
demographic and genetic population parameters. Pre- 
dicted equilibrium population sizes ranged from 30 to 50 
reproductive individuals (average number of reproductive 
plants within G. iaspicula populations = 29.6). The model 
also predicted negative inbreeding coefficients (ranging 
from —0.02 to —0.07), and observed values, based on 
allozyme data, ranged from —0.01 to —0.12 (Hoebee 
et al. 2008). 

Implications for population management 

What are the practical implications for population recovery 
programs when the species concerned has a self-incompati- 
ble mating system? Several authors have demonstrated, 
using interpopulation crosses, a direct link between S allele 
diversity and mate availability which means there is real 
potential for demographic recovery through the addition 
of alleles to genetically depauperate populations (DeMauro 
1993; Willi and Fischer 2005; Pickup and Young 2008). 
However, there has been little attempt to rigorously quan- 
tify how much of a fecundity increase is possible for a 
given number of S alleles, nor has the trade-off with demo- 
graphic approaches (e.g., pollen augmentation, promotion 
of germination, elimination of competition) to increasing 
recruitment and survival been examined. These points are 
illustrated in Figure 6, which shows that when S allele 
diversity is very low, adding even a few new cross-compati- 
ble mating types can significantly increase equilibrium 
population sizes. However, once the number of S alleles is 
above 10-15, purely demographic strategies might be 
better, although the precise point at which one or the 
other approach is likely to be most effective will depend on 
plant life history. Because there is strong negative fre- 
quency dependence with regard to rare mating types in 
self-incompatible systems, new S alleles will be strongly 
self-propagating and thus easy to incorporate into recovery 
programs. 

Summary 

Overall, our results highlight conditions under which 
mate limitation owing to self-incompatibility may be a 
significant constraint, not only on individual fecundity, 
but also for population persistence. However, this is likely 
to be the case for only a small parameter range of life his- 
tories (low b Q and high d) and under conditions where S 
allele diversity is relatively low. At extremes, when popu- 
lation growth rates are either robust or demographically 
very limited, genetic factors probably have only a small 
role to play in determining population viability. Our 
results also show that there can be direct links between 



genetic change and demographic performance in situa- 
tions where self-incompatibility and S allele richness are 
important. Under these conditions, genetic remediation 
can be an effective management approach. 

A crucial unexplored problem in conservation biology 
is to develop a predictive framework for identifying situa- 
tions where genetic intervention is likely to be an effective 
strategy. Such situations may be hard to detect - for 
example, our results indicate that many measures of pop- 
ulation "health" (e.g., F IS , variance in fitness) may show 
little response to changes in population size over a large 
range of values, but will drop off sharply when crucial 
thresholds are reached. Effective management depends on 
the ability to identify populations that superficially appear 
to be stable, but may decline rapidly in the future, given 
current trends (Thrall et al. 1997). 

The model described here incorporates some of the 
demographic and genetic complexities of self-incompati- 
bility in a spatially realistic setting. It is encouraging that 
comparison of model outcomes with data from long-term 
empirical studies of natural systems shows general quali- 
tative and quantitative agreement. However, a number of 
issues warrant further study. First, local spatial genetic 
structure (and its dependence on the extent of pollen flow 
and seed dispersal) will be an important determinant of 
population dynamics and viability (Wells and Young 
2002). Second, our model assumes simple linear domi- 
nance relationships among S alleles, but in nature, these 
interactions are often considerably more complex, poten- 
tially leading to quite different impacts of S allele richness 
on mate availability. Third, most empirical studies of nat- 
ural systems indicate at least some potential for dissolu- 
tion of SI systems (Busch and Schoen 2008), which raises 
the general issue (especially in the gametophytic case) of 
trade-offs between mate limitation and inbreeding depres- 
sion which our analyses do not consider. Finally, our 
results indicate that the occasional addition of S alleles 
could be critical, thus pointing toward the importance of 
metapopulation processes and the potential for "genetic 
rescue" through interpopulation movement of seeds and 
pollen (Thrall et al. 2000). Such genetic rescue effects 
have been demonstrated for inbreeding (Richards 2000) 
and more recently suggested for S alleles (Wagenius et al. 
2007). 
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